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Abstract 

Let k > 0 be an integer and Y a standard Gamma(fc) distributed random variable. 

Let X be an independent positive random variable with a density that is hyperboli- 
cally monotone (HM) of order k. Then Y X and Y/X both have distributions that are 
generalized gamma convolutions (GGCs). This result extends a result of Roynette et 
al. from 2009 who treated the case k = 1 but without use of the HM-concept. Appli¬ 
cations in excursion theory of diffusions and in the theory of exponential functionals 
of Levy processes are mentioned. 
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1 Introduction 

A generalized gamma convolution (GGC) is a limit distribution for sums of independent 
gamma distributed random variables (rvs). The GGCs were introduced by the actuary O. 
Thorin in 1977 when he tried to prove that the lognormal distribution is infinitely divisible 
(see ([22])). He used a technique that later on led the second author of this paper to 
introduce in [6] the concept of hyperbolic complete monotonicity (HCM). The simpler 
concept of hyperbolic monotonicity (HM) was mentioned in (6) pp. 101-102] and more 
carefully studied in [7]. 

The GGCs have got applications in many different fields including infinite divisibility 
(e.g. Steutel and van Harn 12U). mathematical analysis (e.g. Schilling et al. m, stochastic 
processes (e.g. James et al. [T3] and Behme et al. i), and financial mathematics (e.g. 
Barndorff-Nielsen et al. DP)- 
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In 2009 Roynette et al. na proved a novel GGC result that has provided stimulus to 
the present work. In our terminology, they showed that the product of an exponentially 
distributed rv Y and another independent rv X has a GGC distribution provided that the 
density of X is HM. We will give a new and more transparent proof of this result and 
generalize it considerably to cover gamma distributions. 

The paper is organized as follows. In Section [2] the HM, HCM and GGC theory is 
briefly recalled. In Section [3] the main result that the product of a gamma variable with 
shape parameter k and an rv with HMfc distribution has a GGC distribution is given. This 
result can be formulated in several alternative ways. It has also an important extension. 
The proof contains some surprising elements. Applications, analytical as well as stochastic 
process related ones, are given in Section [IJ Finally, in Section 0 some open problems are 
mentioned. 


2 Background 

Basic facts on hyperbolic monotonicity (HM) and generalized gamma convolutions (GGCs) 
are presented here. They are taken from Bondesson pp. Much information about GGCs 
and hyperbolic complete monotonicity (HCM) can also be found in the book by Steutel 
and van Harn m- 

2.1 Hyperbolic monotonicity 

Let / be a nonnegative function on ( 0 ,oo). Consider, for any fixed u > 0, the function 
f(uv)f(u/v), v > 0. Obviously it is invariant under the transformation v K > v~ l . It follows 
that it is a function h(w) of w = v + ir 1 since the value of w determines the set {v, v -1 }. 

Definition 1. A nonnegative function / on (0, oo) is said to be hyperbolically monotone 
(HM or HMi) if, for each fixed u > 0, the function h(w) = f(uv)f(u/v) is non-increasing 
as a function of w = v + v -1 . More generally, it is called hyperbolically monotone of order 
k (HM fc ) if (—iyh^(w) >0, j — 0,1,..., k — 1 and (—1 ) fc-1 /i( fc-1 )(t/;) is non-increasing. If 
this holds for all k > 1, / is also called hyperbolically completely monotone (HCM). 

The class of HM fc -functions is also denoted HM fc . Obviously 

HCM = HM^ C • • • C HM 3 C HM 2 C HM 1 = HM. 

Simple examples of HCM-functions are provided by (with 7 e M, c > 0) : x 7 , e~ cx , and e~ c ^ x . 
It is apparent that the HM^-class is closed with respect to multiplication of functions. For 
/ e HMfc, obviously f{uv)f(u/v) < (/(u)) 2 . It easily follows that log f(e x ) is concave and 
hence that f(x) < Cx 7 for some constants C > 0 and 7 6 1 (depending on /). Every 
HMfc-function / can therefore be modified to an HMfc probability density function (pdf) by 
multiplication by a factor exp (—Six — 5 2 x~ 1 ) (with ^ > 0 and h 2 > 0 arbitrarily small) and 
a normalizing constant. In this paper we are mainly concerned with pdfs. 
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Example 1. Let / be a pdf on (a, b ) C (0, oo) of the form f(x) — C(x — a) a_1 (6 — x)P~ l , 
where C is a constant. It can be shown that for a > 1 and (3 > 1, / is HM& for fc = 
min ([a], [/?]), where [■] denotes integer part. However, if a = 0, then for any value of a, / is 
HMfc for k = [/?]. In particular, the [/(a, 6 ) density is HMi. In fact, in this case it is easy to 
see that h(w) = f{uv)f{u/v) is 0 for all u sufficiently large or small and that for the other 
values of u, h(w ) equals 1 if w is below some bound and otherwise 0 . □ 

Example 2. Let X = U 1 U 2 ■ ■ ■ £4, where the random variables (rvs) Ui are independent 
and uniformly distributed on (0,1). Since — logX has a Gamma(fc, l)-distribution, fx{%) = 
logx) fc_1 ,0 < x < 1. This pdf is HM fc . In fact, h{w) = f{uv)f{u/v) = 0 for 
u > 1 , whereas, for u < 1 , h(w) oc ((logn ) 2 — (logu ) 2 ))^ 1 if u < v < u l (i.e. if 2 < 
w < u + u^ 1 ) and otherwise vanishes. The HM*. result then follows from the fact that 
d(\ogv) 2 /dw is completely monotone (CM). In fact, this derivative can be shown to be 
equal to J 0 °°( 1 + t 2 + tw)~ 1 dt. □ 

The following result, which concerns powers, products and ratios of rvs, is important. Its 
proof (in [Tj) is far from trivial. A main idea in the proof is to use hyperbolic substitutions 
of the form x = uv,y = u/v in certain double integrals. 

Proposition 1. Let X and Y be independent rvs with HM^- densities (X ~ HM*,,!' ~ 
HMfcj. Then, for any q G M with \q\ > 1, we have X q ~ HM^. Moreover, X-Y ~ HM^ and 
X/Y ~ HM fc . 

A simple consequence of Proposition Q] (with one of the rvs exponentially distributed) 
is that the Laplace transform of an 11M/. function is 11M k . Let A" have the HM 2 -density 
f[x) = 2max(0,1 — x) and let E~f/(0,1) (with an HMi-density). Then it can be shown 
that X/Y rfj HM 2 . Thus there is no trivial extension of Proposition [fl 

The HMi-densities can be identified as follows (see 0 )- 

Proposition 2. We have X ~ HMi if and only if Y = log A" has a pdf that is log- 
concave, i.e. log fy(y) is concave. Equivalently, X ~ HMi if and only if fx{x) = 
C exp(— dy ), where if is non-decreasing, C a constant, and x 0 is suitable chosen. 

With this, the well-known fact that logconcavity is preserved under convolution (see 
e.g. m pp. 17-23]) becomes a simple consequence of Proposition Q] for k — 1. 

Typical HCM (=HMoo) pdfs have the form f(x) = Cx /3_1 n”=i(l + c^) _7 % where the 
parameters are positive, or are limits of such densities. In fact, all HCM-densities (and 
functions) are such limits. An open problem is to find canonical representations for HM^- 
densities for 1 < k < 00 . 

The HMfc-class of densities (functions) can alternatively be described by the condition 
that 

h{w) = f(uv)f(u/v) = c u + f (A - w) k ~ l H u (d\), (1) 

J (lO,Oo) 

where c u > 0 and H u {dX) is a nonnegative measure. The simple example f(x) = x 7 gives 
c u = u 2,1 and H u (d\) = 0. However, for a pdf we must have c u = 0. The representation (JT]) 
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follows from a representation of the non-increasing function (—1 ) k ~ 1 h^ k ~ 1 'l (w) as an integral 
over (w, oo) (or possibly [w, oo)) of a nonnegative measure. For instance, for k = 2 we put 
— h'(w ) = f 1 (w < X)H u (dX). We then get, by a change of the order of integration, 

l(tc < w < X)dwH u (dX) — I (A — w)H u (dX). 

J (w, oo) 

The representation dTJ) was derived and used in [7j. For functions with monotone deriva¬ 
tives up to some order it seems to have been first used by Williamson [23] . 


h(w) — h(oo) = — h'(w)dw = 


2.2 Generalized gamma convolutions 

Convolving different gamma distributions, Gamma(ii,l), with pdfs and Laplace transforms 
(LTs) of the forms f(x ) = (T(u))~ 1 x u ^ 1 t u exp(-xt) and 0(s) = respectively, and 

then taking weak limits, Thorin [22] was led to the following definition. 

Definition 2. A generalized gamma convolution (GGC) is a probability distribution on 
[0, oo) with LT of the form 

0 (s) = exp ^-as + ^ log U ( dt )j > 

where (the left-extremity) a > 0 and U(dt ) is a nonnegative measure on (0, oo) (with finite 
mass for any compact subset of (0, oo)) such that ^ | log t\U(dt) < oo and f (l ^ t~ 1 U(dt) < 
oo. 

The GGC-class of distributions is closed with respect to (wrt) addition of independent 
random variables and wrt weak limits. Each GGC is infinitely divisible and each convolution 
root of a GGC is a GGC as well. The pdf f(x) of a GGC is strictly positive on (a, oo) 
and, if a = 0 and /3 = f (Q ^ U(dt) is finite, then f(x) = x^^hl^x), where h(x) is completely 
monotone (see 0 p. 49]). 

The pdf of a GGC need not be HMj. For instance, for a gamma distribution with shape 
parameter less than 1 and shifted to have left-extremity a > 0 the pdf is not HM 4 . An 
HMfc-density, which may have compact support, is in general not a GGC. However (see 0 
Theorem 5.1.2]): 

Proposition 3. If the pdf f on (0, 00 ) is HCM, then it is a GGC. Thus HCM C GGC. 

Many well-known pdfs are HCM and therefore also GGCs and hence infinitely divisible. 
For instance gamma densities are HCM. Then it follows from Proposition [T] (for k = 00 ) 
that also the power q, q > 1, of the ratio of two independent gamma variables has a density 
that is HCM. This density is of the form f(x) = Cx^~ 1 (l + cx a )~' y ,x > 0, with a = q~ l . 
Every lognormal density is also HCM. 

The next proposition gives a characterization of the LT of a GGC (0 Theorem 6.1.1]). 

Proposition 4. A function <p(s) on (0, 00 ) is the LT of a GGC if and only if 0(0+) = 1 
and 0 is HCM. 
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This result will be our basic tool in Section 4. Since the LT of an HM(. function is HM*,, 
and this also holds for k = oo, Proposition [3] can be seen as a consequence of Proposition [4j 
Using another complex characterization of the LT of a GGC, we can get the following result 
(IS Theorem 4.2.1]). 

Proposition 5. Let Y ~ Gamma(l, 1) and let X > 0 be an independent rv with a density 
f(x) that is logconcave (or only such that xf(x) is logconcave). Then Y/X ~ GGC. 

One should notice that in Proposition [5] the rv X is not assumed to have an HMi-density. 

Proposition 6. If f(x) is the density of a GGC and x~ a f(x), where a > 0, can be nor¬ 
malized to be the pdf g(x) of a probability distribution, then g{x) is also the pdf of a GGC. 

This result is only a limit case of [6) Theorem 6.2.4]. The following recent result from 
[S] needs to be mentioned. It can be proved by the help of Proposition [4J 

Proposition 7. Let X ~ GGC and Y ~ GGC be independent rvs. Then X -Y ~ GGC. 

Well-known examples of GGC distributions include the log-normal distribution and 
positive strictly a-stable distributions. Also, each negative power of a gamma variable is 
shown to have a GGC-distribution in [!9]. Bosch and Simon ra and Jedidi and Simon [14] 
give other novel results on HM, HCM, and GGC distributions. 

3 Main result 

Here the main result is presented as a theorem in Section 13.11 Moreover comments are 
given. The proof is presented in Section 1X21 

3.1 Formulation of the main result and comments 

Theorem 1 . Let k > 1 be an integer. Let Y ~ Gamma(fc, 1) and X ~ HM& be independent 
rvs. Then Y-X ~ GGC and Y/X ~ GGC. 

We give some comments on the above theorem. 

Remarks 1. (i) For k — 1 Theorem |T] differs from Proposition [5j One should notice that 

X ~ HMi <;=> X~ l ~ HMi but logconcavity of fx is not equivalent to logconcavity 
of fx- 1 - One can also notice that every gamma density is HCM (and thus HMi) but 
only logconcave when the shape parameter is > 1. 

(ii) In the case k — 1 the LT <f>i (s) = f^°(x + s)~ 1 xfx{x)dx of Y/X for independent 
Y ~ Gamma(l, 1) and A" ~ HMx is the Stieltjes transform (or double Laplace trans¬ 
form) of the measure xfx(x)dx. For k > 1 the LT 0 fc (s) = J 0 °° (x + s)~ k x k fx(x)dx 
coincides with the so-called generalized Stieltjes transform (of order k ) of the measure 
x k fx(x)dx. In that sense the above theorem can be restated as follows: Assume fx(x) 
is an HM*, function. Then the k-th order generalized Stieltjes transform of x k fx{x)dx 
is HCM, i.e. it is the LT of a GGC. 
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(iii) Clearly Theorem [Tj remains true if Y ~ Gamma(fc, 6) for any 9 > 0, since in this 
case 9Y ~ Gamma(fc, 1). Considering Y ~ Gamma(/c, k) and letting k —> oo we get 
that Y — y 1 in probability. Hence for X ~ HM& with k fixed it is neccessary in the 
theorem to have a restriction upwards on the shape parameter of Y since otherwise 
it would incorrectly follow that HM*. C GGC. For instance, if Y ~ Gamma(2,1) and 
X ~ 17(1,2), then f x is HMj but Y/X / GGC. 

(iv) Letting again k —> oo and so that Y —> 1 in probability, we get back Proposition [3] as 
a limit case of Theorem [TJ Since a Gamma(fc, 1) density is HCM, it also follows that 
the class of GGCs provided by Theorem [T] is closed wrt multiplication and division of 
independent rvs. However, if Z = Y-X with Y ~ Gamma(fc, 1) and X ~ HM&, it is 
not true that Z _1 always has the same representation. 

(v) Theorem[T]can also be expressed in the following way. Any scale mixture of Gamma(fc) 
distributions with a scale mixing HM^-density is a GGC. It is well known m 
Theorem 3.3, p. 334]) that any scale mixture of Gamma(l) distributions is infinitely 
divisible (ID). More generally, any scale mixture of Gamma(2) distributions is ID 
m)- However, for k > 2 ID fails to hold in general for such mixtures. 

There is a nice extension of Theorem Q] which we see as a corollary of it. 

Corollary 1. Let Y ~ Gamma(r, 1) be independent of X ~ HM& where r > 0 and k is an 
integer such that k > r. Then Y-X ~ GGC and Y/X ~ GGC. 

Proof. Since X ~ HM^ if and only if 1/X ~ HM&, it suffices to consider the ratio Z = Y/X. 
Let a = k — r and let Y' ~ Gamma(fc, 1). Then 

POO -t POO 

fz(z)= xf Y (zx)f x (x)dx =—— x{zx) r - 1 e~ zx f x {x)dx 
Jo 1 \ r ) Jo 

= ^rr\z~ a [ xfyr(zx)x~ a fx{x)dx. 

r ( r ) Jo 

Here x~ a fx{x ) is HM^ and so is, for any 5 > 0, x~ a e~ 5 / x fx{x). Letting if necessary 5 > 0 
and normalizing this latter function to become the pdf of an rv X', we get from Theorem Q] 
that Y'/X' ~ GGC. Using then Proposition [ 6 ] and letting S —* 0, we conclude that 
Y/X - GGC. □ 

3.2 Proof of the main result 

The proof of Theorem |T| is given in two parts. First the case k — 1 is treated. This 
proof is short but contains the essential ideas. The proof in the general case becomes more 
technical. Of course, we use the HCM-characterization of the LT of a GGC and hyperbolic 
substitutions in the proofs. For the transformation T — t + £ -1 , we avoid to use the inverse 
transformation t = T/2 ± a/T 2 /4 — 1. In fact, the HCM-concept was introduced in the 
early 1990s in order to avoid, at least in presentations, such inverse transformations. 
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Proof of Theorem [Jl k—1. It suffices to consider the ratio Y/X, where Y ~ Gamma(l, 1). 
The LT (f)(s) of the distribution of the ratio is given by, with / = fx, 


<p(s) = U(exp(— sY/X)) = / E(exp(—sY/x))f(x)dx = 

Jo 

For fixed s > 0, consider 


x 


X + s 


f(x)dx. 


POO poo 


J = 4>(st)4>(s/t) = 


xy 


I o Jo (x + st)(y+ s/t) 


f(x)f(y)dxdy. 


In view of Proposition [IJ we only have to show that J is completely monotone (CM) wrt 
T = t + £ -1 . We make the hyperbolic substitution x = uv,y = u/v with Jacobian with 
modulus 2 u/v. Hence 


J 



2 u u 2 

v (uv + st)(u/v + s/t) 


f(uv)f(u/v)dudv. 


Using the representation f{uv)f{u/v) = J^ woo ^H u (dX), where H u (d\) is a nonnegative 
measure and w — v + v ~ l , letting b = 6 (A) > 1 be such that 6 + b~ l = A, letting a = u/s, 
and changing the order of integration, we get by some simple algebra that 


J 



Ii/b (v + t/a){y + at) 


dv ) H u (dX)d'u. 


=.Ji 


It is now evident that it suffices to show that for each 6 > 1 and each a > 0 the interior 
u-integral .J\ is CM wrt to T = t + £ -1 . For 6 = 1, J\ — 0, so it suffices to consider the case 
6 > 1. The integral J\ is a function of T since the change t n- £ _1 leaves J\ invariant which 
is shown by the substitution v = 1/V. Now J\ can be calculated explicitly. In fact, by a 
partial fraction expansion we have for a ^ 1, 


t 

(v + t/a)(v + at) 


a — a 1 \v + t/a v + at 


and hence, for a ^ 1 , by an integration and some simplification, 


1 f [t + ab){t + ( 06 ) X )A 

— ° g \(t+ a/b)(t+ 6 /a) J 



log 


T + A 
T + B ’ 


where A = ab + (a 6 ) _1 ,H = a /6 + 6 /a. For a = 1, J\ — (6 — 6 ^ 1 )/(T + 6 + 6 ^ 1 ). Since 
a n- a -1 leaves J\ unchanged, we may without restriction assume that a > 1 (and as earlier 
6 > 1), and then A > B and J\ > 0. Moreover, we get that the k- th derivative of Ji, i.e. 
here its first derivative, has the form 


dJ\ 1/1 1A 1 B - A 

~dT = a-a- 1 \T + A ~ T + B J = a - a " 1 {T + A)(T + B) 
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and this derivative is negative. Since (T + A) l (T + B) 1 is CM, we get as desired that 
(—1 yj[ j \T) > 0, j = 0,1,2,..., and the proof is complete. □ 

We now proceed with the general proof of Theorem |T] for any integer k > 1. We shall 
see that the above proof needs some complementary arguments. 

Proof of Theorem [3 general k. Let Y ~ Gamma(fc, 1) and X ~ HM& be independent. 
Then the LT f>(s) of the distribution of Y/X is given by 


0(s) = 


x 


X + s 


f(x)dx. 


Hence using (|T|) 


J = <f(st)<f(s/t ) = 


2 u 


h k 

x y 


u 


(x + st) k (y + s/t) k 

2k 


f (x) f (y)dxdy 


v (uv + st) k (u/V + s/t) k 


f(uv)f(u/v)dudv 


2 u 


u 


2k 


(A — w) k 1 H u (d\)dudv , 


'0 

,-i 


v (uv + st) k (u/v + s/t) k 


where w = v + r; _1 and H u (d \) is a nonnegative measure. Again we let b = 6(A) > 1 be 
such that b + 5 -1 = A and put a = u/s. After a change of the order of integration with the 
v —integral as the inner integral and noticing that b + b~ x — v — u _1 = (b — v)(v — b~ l )/v, 
we see by some algebraic manipulations that it suffices to show that the integral 


-i\\fc-i 


J k = hdv , where I k = 
Ji/b 


t k ((b — v)(v — b x )) 
(u + t/a) k (v + atf 


( 2 ) 


is CM wrt T — t + t~ x . The same argumentation as in the case k = 1 shows that ,J k is a 
function of T. An important fact is that J k can be calculated explicitly for all integers k 
although J/,. becomes complicated for k large. Since I k is a rational function of v, we can 
get an expression for J k by using first a partial fraction expansion of I k wrt v. However, it 
is more efficient to use an alternating generating function: 


OO oo 

GF(z) = = / E(- 

fc=i Jl /b k=1 


z) k ~ l I k dv 


1 


-dv. 


Ji/b ( v + t/a)(v + at)/t + z(b - v)(v - 6 _1 ) 

Minimizing over v and t, we see that the series is absolutely convergent at least if \z\ < 
Since the denominator in the integrand is a quadratic function of v and as such 


(l+a) 2 b 

a function can be factorized into two real linear factors for z > 0, we get by integration and 
considerable simplification with the notation a = a + a -1 , f3 = b + 6 _1 that 

1 


VA 


GF(z) 


log R, where 










A = (a + pz ) 2 - 4 - Az 2 + 4^T and i? = 


T — 2z + |/5(a + /fo) + |(6 — b 1 )-\/A 
T — 2z + \P(a + Pz) — \(b — b~ l )\fK 


It is far from obvious but some calculation shows that 

d, o 2(6-6- 1 ) 

dz l0g ' 

In fact, the product of the numerator and the denominator in R does not depend on z so 
the derivative above is just twice the derivative of the logarithm of the numerator in R. Of 
course, Jk = Jk(T ) = GF^ k ~ l \ 0). Now it is not difficult to see that with, as earlier, 

A = ab + ( ab ) _1 and B = a/b + b/a, we get 


J k (T) = P k (T) + Q k (T) log 


T + A 
T + B 


(3) 


where P k (T) and Q k (T) are polynomials in T of degrees k — 2 and k — 1, respectively. For 
k — 1, P k (T) vanishes. For k — 1,2, and 3, we have 

Pi(T) = 0, P 2 (T) = - P 3 (r) = -3 5 y|y I (2r + /i + B), 

Qj(T) = -Vr. Q 2 (T) = g±4t# , Q 3 (T) = 

By using the above expressions for P k (T ) and Q k (T) one can easily verify that at least for 
k = 1, 2, 3 we have somewhat surprisingly 


d% 

dT k 


= (-l)\k-l)\ 


(i b - b~A 


-l\2k-l 


(T + A) k (T + B) k ' 


(4) 


To see that (J4|) is completely general, some additional argumentation is needed. Since P k (T ) 
has degree k — 2, it has no influence at all on the fc-th derivative of J k . Since Q k (T) has 
degree k — 1 and hence Q^\T) = 0, it also follows from (J3J) after some reflection that 


Rk(T) 


d k J k _ 

dT k (T + A) k (T + B) k ’ 


(5) 


where Rk{T) is a polynomial of degree at most 2k —1. To see that really Rk{T) is a constant, 
(— l) k (k — 1 )!(6 «» 6 ~ 1 ) 2fe_1 , we look at the case when t —y oo. Then T = t + t^ 1 is very close 
to t. From (j2J) we get that 


J k ~ ^ j^((b-v)(v -b x )) 


_i\\ k —i 


(6 — 5-i)2fe-i 

du ~ B(k,k) -—- as too, 

T k 


where B(k , /c) = (A; — l)!(fc — l)!/(2/c — 1)!. Since k(k + 1) • • • (2k — l)B(k, k) = (k — 1)!, we 
also get that 

.( b - b - 1 ) 2 *- 1 


d fc A, 


dT fc 






as t —> oo. 
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Since Rk(T) in (151) is a polynomial, this asymptotic relation can only hold when the poly¬ 
nomial is a constant and hence (J3J) holds for all k. To complete the proof, we use that 
( T + A)~ k {T + B)~ k is CM and hence (— 1 Y J^\t) > 0 for j = k, k + 1,.... Then it only 
remains to verify that these inequalities also hold for j — 0 , 1 ,..., k — 1 . Using the same 
argumentation as above we have that, for each j > 0, J^\t) = 0{T ~ k ~as T —* oo. In 
particular J^\T) vanishes at T — oo. It follows that 

/ OO 

4 J+1) (f) df, j =0,1,2,... 

We see that the sign of is opposite to that of j[ k \T). The same then holds for 

the sign of jj^ 2 \t) compared with that of jjf ^(T), etc. This shows that Jfc(T) is CM 
wrt T as desired. □ 

Remark 2. In the general case some technical details have been omitted in the proof, 
ffowever, it is easy to check all statements by using a program for symbolic algebra. In 
fact, the simple form in ([4]) for the derivative J^(T) was discovered in that way. 


4 Applications 

Theorem Q] and Corollary [Tj have a wide range of possible applications. We will discuss a 
few in this section. 


4.1 Excursion theory 

The random process foundations for the research carried out in this article have been laid 
by Roynette et al. HU and Salminen et al. [T8] . In these articles the authors study excursion 
times of recurrent linear diffusions on M + . More precisely, given an M + -valued recurrent 
diffusion (X t ) t > 0 and defining the last and the next visit in 0 via 


9t := sup{s < t,X s = 0}, 


d t := inf{s > t, X s = 0} 


they are interested in the rvs 


y(i) 

v 




y( 3 ) 

p 


dz p — 9z p , 


( 6 ) 


where Z p denotes an exponential rv with density pe pz , z > 0, independent of (X t ) t > 0 . In 
[18j it is shown, that all Y W are infinitely divisible, while in Id the authors give conditions 
for yW to have GGC distributions. These conditions are stated in terms of the Krein 
measure of the Levy measure of the inverse local time at 0 of (X t ) t > 0 . 

For their proof of the GGC property, Roynette et al. first show, for k — 1, a reformulation 
of Theorem mam Theorem 2]). They do not use the HM-concept but define a class C 
of functions which essentially coincides with the class HMi. The proof of DU Theorem 2] 
then relies on the HCM-characterization of the LT of a GGC (Proposition [4]). Although 
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also our proof for k — 1 uses Proposition 0 it is shorter than theirs because of our use of 
a suitable hyperbolic substitution in a double integral and the avoidance of certain inverse 
transformations. 

Further in [1TJ the LTs of the Y^s are shown to be Stieltjes transforms of measures 
whose densities are LUvR (compare with Remark din)). 

In the following we will indicate via an example how one can also use our main theorem 
in the case k — 2 to prove the GGC property of Yp 3 ^ as defined in (151) . Therefore we briefly 
recall some notation from [18] and HU- 

Let 0 be the continuous local-time of {X t ) t > 0 at 0 and (t u ) u > 0 its right-continuous 
inverse. Then ( r u ) u > 0 is a subordinator and as such has a Levy exponent -0 and a Levy 
density u, i.e. 

/ /'OC 

E[e~ XTu ] = e~ u ^ x) = exp ( ~-u J (1 - e~ Xx )u(x)dx 
where further v has the Krein representation 

poo 

u(x) = / e~ xz K[dz ) 

Jo 

with Krein measure K of u. 


Proposition 8. Assume that the Krein measure K is such that the function f(u ) defined 
via 

f(u~ 3 ) = / K(dz) 

J (u —p) VO 

is an HM 2 function. Then Yp 3) ~ GGC. 


/o\ 

Proof. It was shown in [18], that the distribution of Yf is a Gamma(2)-mixture. In 
particular it can be deduced from [18], Equations (46), (49) and (50)] that the density of 
Yp 3> is given by 


/ v (3)(u) — 


if(p) 


ue 


Ii(dz)dx 


r x—p 


(o\ 

which shows that Yf = Y ■ X where Y ~ Gamma(2,1) and X is independent of Y with 
density fx dehned via fxiu" 1 ) = f™_ p K(dz). Thus the claim follows from Theorem Q] 
in the case k = 2. □ 


4.2 Exponential functionals of Levy processes 


Let £ = (£ f )t >0 be a Levy process such that —> —00 as t —> 00 . Then the exponential 

functional of f is dehned as 


T := / e^dt. 

J (0,oo) 


Such exponential functionals appear as stationary distributions of generalized Ornstein- 
Uhlenbeck processes and they have attracted a lot of interest throughout the last years (see 
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e-g. HQ, the survey paper p] or the more recent contributions PI El US] to name just very 
few references). 

ft is known, that Jg ~ GGC in several cases. E.g. Dufresne (e.g. [5], Equation (16)]) 
showed that Ij = -^G^a/a 2 w h ere G 7 ~ Ganimacy, 1), whenever £ is a Brownian motion 
with variance a 2 and drift a < 0. Concerning processes £ with jumps, one has for example 
the following Proposition. 

Proposition 9. Suppose that £ is a compound Poisson process, i.e. £ t = 'Ylnl l X i with 
i.i.d. jump heights X i: i = 1,2 ,..such that —oo < E[Xf\ < 0 and e Xl ~ GGC. Then 
- GGC. 

Proof. The proof can be carried out along the lines of the proof of [3> Proposition 3.2] using 
the more recent Proposition [7] □ 


Still, assuming that £ t = at—N t , t > 0 , for a < 0 and a subordinator (i.e. a nondecreasing 
Levy process) (N t )t>o, one easily observes that It has bounded support and therefore cannot 
be infinitely divisible such that in particular Jg ^ GGC. 

In [16] , based on the Wiener-Hopf factorization of Levy processes, the authors obtain 
factorizations of exponential functionals. In particular, in case of a spectrally negative 
process £ with —>• — oo, they prove that 


k 


d 


Gf/ 


where H = (H t )t>o is the descending ladder height process of £ and G 7 ~ Gamma(y, 1 ), 
with 7 depending on the characteristics of £, is independent of H. We refer to [T9] or |Tj 
for any further information on Levy processes, their characteristics and their Wiener-Hopf 
factorizations. 

Since —H = (—H t )t>o is a subordinator with drift an and Levy jump measure uh, say, 
it follows from m Example B] that if H is noil-trivial, then In admits a density f ( 5 ) which 
fulfills the integro-differential equation 


(1 - a H s)f(s ) = / u H (log(t/s))f(t)dt, 


where u H (x) = Vh((x, oo)). 
authors prove that 


In particular, if v H (s) = ce bs , b, c > 0, and an > 0, the 



Zb+l,c/a,H i 

a H 


where Z a ^ is a Beta rv on (0,1) with parameters a,/3 > 0. Hence Ih ~ HM^ for k < 
min ([6 + 1 ], [ c/a H ]). 

Now by Corollary |TJ I^ is the reciprocal of a GGC if k > 7 . Notice that in general 
such inverses of GGCs are not GGCs themselves. However, in this case L ~ 11M 1 .. since 
Gl ~HCM. 

Conversely, if again uh(s) = ce ~ bs , b, c > 0, but an = 0, then Ih itself is Gamma 
distributed and so I g ~ GGC. 
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4.3 Constructing GGCs 

Using Theorem |Tj we can construct explicit densities and LTs of GGCs as we shall do in 
the following. 

Examples 3. (i) Let Y ~ Gamma(l, 1) and X — U ~ U(0,1) be independent. Then 

we have the following LTs and pdfs for YU and Y/U, respectively: 

<hu{s) = lQb ^ + — , <t>Y/u{s) = 1 + s log(-^—), 

s 1 + s 

POO 1 

fvu(x) = Ei(x ) = / y~ l e~ yx dy, fv/u(x) = —(1 - (1 + x)e~ x ). 

J i x 

By Example [2] we have U ~ 11M, and hence by Theorem |T] the above pdfs are GGCs 
and the LTs are HCM. 

(ii) Now let Y ~ Gamma(2,1) and X = min(f/i, U 2 ), with U\, U 2 ~ U( 0,1) independent 
and independent of Y. Then fx(x) = 2(1 — x),0 < x < 1, which belongs to HM 2 . 
We can also represent X as X = UiU^ 2 . We get the following LTs and pdfs: 

(Pyx(s) = -(1 - 1 ° g ^ 1 + g ^ ), 4>y/x(s) = 1 + 6s + (6s 2 + 4s) log(-^—), 

s s 1 + s 

fvx(x) = 2e~ x - 2 xEi(x), fv/xix) = -\(-12 + 4x + (2x 2 + 8x + 12)e~ x ). 

x 6 

Again by Theorem [T| the pdfs are GGCs and the LTs are HCM. 

Many similar examples can be obtained from Corollary Q] by letting Y ~ Gamma(r, 1) 
with a real r. 


5 Final comments 

There are reasons to believe that Theorem [Tj (as well as Proposition [I]) can be extended to 
cover the case that k is any real number > 1. Maybe it can even be extended to all real 
k > 0. As a definition of an HM^ function in the real case, the integral representation (JTJ) 
can be used. For any real j and k such that 0 < j < k we have HM*. C HMj. To see this, 
one can use (JTJ) together with the simple formula 

For k > 1 the HM*. class is closed wrt multiplication of functions. However, it is not closed 
for k < 1 which the example f(x) = (1 — x )^ 1 / 2 illustrates. The technique which we have 
used to prove Theorem |T] for integers k cannot be applied in the general real case since 
it much depends on an explicit calculation of the integral Jk in ([2]). However, numerical 
experiments indicate that Jk is CM as a function of T = f+ f _1 for all k > 0. An important 
problem for the future is to prove that so is the case. 
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Let A and B denote classes of probability distributions. We denote by A x B the class 
of distributions generated by Y ■ X for Y ~ A and A" ~ B with Y and X independent. 
Theorem Q] and Proposition [7] can then be formulated as Gamma(fc) x 11 M p. C GGC and 
GGC x GGC C GGC, respectively. 

One may wonder about the largest class Tik such that Gamma(fc) x TLk C GGC. Appar¬ 
ently, because of Theorem |Tj and Proposition [3 77/. D 11M/, x GGC. One may also wonder 
about the largest class Qk such that Qk x IIM/,. C GGC. Of course, Qk D Gamma ( A: ) x GGC. 
Could possibly Qk D GGC(fc), where GGC(A:) denotes all GGCs with left-extremity 0 
and total [/-measure at most k? To prove this possible result, it suffices to show that 
Y ■ X ~ GGC when X ~ I IM/. and Y is a finite sum of independent gamma variables with 
a shape parameter sum not exceeding k and varying scale parameters. 
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